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Abstract 

In this paper we use our recent technique for estimating the turbulent component 
of the magnetic field to derive the structure functions of the unpolarized emission as 
well as that of the Stokes Q and U parameters of the polarized emission. The solutions 
for the structure functions to 350-/mi SHARP polarization data of OMC-1 allow the 
determination of the corresponding turbulent correlation length scales. The estimated 
values for these length scales are 9?4±(X'l, 7 / . / 3±0 / /l, 12?6±(X'2 (or 20.5±0.2, 16.0±0.2, 
and 27.5 ± 0.4 mpc at 450 pc, the adopted distance for OMC-1) for the Stokes Q and 
U parameters, and for the unpolarized emission N, respectively. Our current results 
for Q and U are consistent with previous results obtained through other methods, 
and may indicate presence of anisotropy in magnetized turbulence. We infer a weak 
coupling between the dust component responsible for the unpolarized emission N and 
the magnetic field B from the significant difference between their turbulent correlation 
length scales. 

Subject headings: ISM: clouds - ISM: magnetic fields - polarization - turbulence 

1. Introduction 



Magnetic fields play a crucial role in the formation of stars through var ious processes, r ang- 
ing from magnetic support against gravitational collapse (see the reviews by I Shu et al .11 19991 . and 
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Mouschovias Sz Ciolekl Il999l ) to magnetic braking (|Nakanolll984l ). On the other hand, turbulence 
has also been sugge sted to be a regulating fac tor i n the development of the star formation process 
(see the reviews by iMac Low fc Klessenl |2004| . and lElmegreen fe Scald 12004 ) . The re lative impor- 
tance of these two agents d uring the star formation process remains controversial (e.g.. iMouschovias 
20091 : ICrutcher et al.ll201oh . 



The magnetic field threading a molecular cloud is composed of a large-scale ordered component 
and a turbulent (or random) one. Cha racterization of the latter, e.g., through dust po larization 
measurements, is highly desirable (e.g., lLai et al.ll200ll , 2002, 2003 . ICrutcher et al.l I2Q04T) . For ex - 
ample, with this information one can use the method introduced bv lChandrasekhar fe Fermil (|l953l ) 
together with the estimates of gas density and line of sig ht velocity dispersion in mol ecular clouds 



(from appropriate molecular species, e.g., H CO + ; see iHildebrand et al.l I20091 . and iHoude et al. 



2009) to derive the strength of the pla n e of the sky component of the magnetic field in thes e regions 



(e.g., iHoudd 12004 : iGirart et al.l 120061 : ICurran fc Chrvsostomoul 120071 : lAttard et al.1 l2009h . It can 



also provide us with a measure of the turbulent energy content of the gas. 

A potential problem in using such a method is the necessity to make assumptions on the 
structure of the large-scale component of the field. These assumed models for the morphology of 
the magnetic field, even though derived from a polarization map, may still lead to imprecise values 
for the tur bulent component as th e assumptions regarding the large-scale will not follow its true 
structure. Hildebrand et al. (|2009l ) introduced a new technique where the dispersion function of 



the polarization angle is used to estimate the turbulent component of the field without making any 
assumptions on the morphology of the large-scale magnetic field. This method has been studied in 
more details by IHoude et al.l (|2009l ) , where the effect of signal integration was taken into account 
in the evaluation of the ratio of the turbulent to large-scale magnetic field strength. This has also 
allowed, among other things, the determination of the turbulent magnetic field correlation length 
scale. 



In this paper, we will use the methods introduced in IHoude et al.l (J2009J) to compare the tur- 
bulent character of the polarized emission (emanating from the dust component that is coupled to 
the magnetic field) and that of the unpolarized emission by estimating their turbulent correlation 
length scales. We will start by deriving the cloud- and beam-integrated structure functions for the 
Stokes Q and U parameters and the unpolarized emission ./V in Section 2. We will then proceed 
with presenting the solution for the structure functions assuming Gaussian profiles for the beam 
and the turbulent autocorrelation functions of Q, U, and N. In Section 3, we use the previously 
published 350- /xm SHARP polarization map of OMC-1 of IVaillancourt et al.l (120081 ) to determine 
the turbulent correlation lengths, and then use these results to compare the turbulent character 
of the polarized emission and the unpolarized emission and gain information on the nature of the 
turbulence and magnetic field existing in the region of interest. 
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2. Analysis 

2.1. The Cloud- and Beam-integrated Structure Function of the Polarized 

Components of the Emission 

The linear polarized emission is observationally determined through the measurements of the 
cloud- and beam-integrated Stokes parameters 



Q(r) 



H(r - a) 



Q (a, z) dz 



d 2 a 



(1) 



U(v) = J J H(r - a) ±Ju(z,z)dz 



with 



d a, 



(2) 



P(r) = jQ 2 (r)+U 2 (r) 



(3) 



the integrated polarize d emission, H (r) t he beam profile, and A the maximum depth of the cloud 
along any line of sight (jHoude et al.ll2009l ). The unpolarized emission N can be estimated through 
the measurements of the total emission / and the polarized emission P with 



N(r) = I (r) - P (r) 
In analogy to Q and U, N is defined as 



(4) 



N(t) 



H(r 



A 



N (a, z) dz 



d 2 a. 



(5) 



The quantity for N (r) defined in Equation ([5]) is, however, not directly measurable. We therefore 
use the following approximation 



iV(r)~ J(r)-P(r). (6) 

The two dimensional integrals in Equations (HJ. (|2|), and ([5]) are over all space. The z-axis is along 
the line of sight (unit basis vector e^), and r is the two dimensional polar radius vector on the plane 
of the sky (unit basis vector e r ) such that 



x = re r + ze z . 



(7) 
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Due to the mathematical similarity between the definitions of Q, U, and N, the analyses will be 
shown for only one of them (the Q component), the results for U and N follow in a straightforward 
manner. 



Similarly to th e characterization of the magnetic field presented in lHildebrand et al.l (|2009l ) and 



Houde et al.1 (120091 ) . we will divide the emission - either polarized or unpolarized - into ordered and 



turbulent components 



Q (a, z) = Qq (a, z) + Q t (a, z) . (8) 



In order to gain a quantitative estimate of the turbulent component of the magnetic field in molecular 
clouds, it is necessary to make some assumptions about the Stokes parameters (Q and U) as well 
as the unpolarized emission N. More precisely, we assume statistical independence between the 
ordered and turbulent components. We will therefore have the following averages for any two points 
x and y 



(Qo(x)) = Q (x) 

(Q*(x)) = o 

(Qo(x)-Q*(y)> = (Q (x))-(Q t (y)) = 0. (9) 

Moreover, we assume that Q is stationary and isotropic (see Equations |15] and [16] below). Note 
that the distance between the two points where the Stokes parameter Q is measured is confined to 
the plane of the sky unless otherwise stated. 

The autocorrelation function of Q can be introduced as 

R Q (v,u) = (Q( a ,z)Q( a ',z')) (10) 

with u = \z' — z\ and v = |a' — a|. Due to the assumed statistical independence of the ordered and 
turbulent components, it is possible to decompose these parts of the autocorrelation function in the 
following way 

R Q (v, u) = R Qfl (v, u) + R Qit {v, u) , (11) 

with 



R Q,j K u ) = (Qj ( a > z ) Qj ( a '> z ') ) 
in which j stands for "0" or "t" for the ordered and turbulent parts, respectively. 



(12) 
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In order to estimate the structure function of the Stokes parameter, it is necessary to specify 
some characteristic s of the autocorrelation function, as well as the telescope beam profile. Following 



Houde et al.l f)2009l ) we write 



Rq (v,u) = R Q>0 (v,u) + (Q 2 ) e -(- 2 +« 2 )/^ (13) 

with 5q the correlation length scale for the turbulent component of Q. For the time being, we 
assume that U and N have turbulent correlation length scales 8u and 6n, which are potentially 
different from Sq. Moreover, these correlation lengths are both assumed to be much smaller than 
the thickness of the cloud (A). The beam profile is also assumed Gaussian 



H(r) 



1 

2ttW 2 ' 



-r 2 /2W 2 



(14) 



with W the beam radius. 



For the structure function (IFalceta-Goncalves et al.ll2008l : iHildebrand et al.ll2009l ) of Q we start 
with the following definition 



AQ Z (£) 



\Q(t)-Q(t + £) 



Q 2 )-(QQ(t)) 



(15) 



where we have used the aforementioned assumptions of isotropy and stationarity and 



[QQ(t)) = {Q(r)Q(r + £)). 



(16) 



Under these constraints, decomposing the ordered and turbulent compone nts of Q, and incorporat- 
ing the analytical solution provided in Appendix A of lHoude et al.l (|2009l ) we can write 



5 2 Q + 2W 2 ) A 



(17) 



where 



(')> = 2 



H (a) H (a' + 1) 



11 

1_ x) R Qfi ( v ' u ) du 



j2 I j2 

a a a a. 



(18) 



Inserting Equations (I17p and fjl8[) into Equation (1 1 5 [) we get 
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AQ (£) 



[(no)) -(am 



(19) 



+2V2tt (Qf 



5 2 + 2W 2 ) A 



1 -e 



-i 2 /2(5 2 Q +2W 2 ) 



The first term on the right-hand side of the above equation (within brackets) is solely due to the 
ordered component of the Stokes Q parameter and not turbulence. T his term can be ex panded as 
a Taylor series (similarly to the analysis presented in the Appendix of iHoude et al.l 12009) 



«(<>)> -<*(*)>=!> 



;2i 



(20) 



8=1 



while the isotropy in £ is incorporated by performing the summation only on terms of even orders 
in £. For small displacements satisfying I < W, the ordered term can be described adequately by 
keeping the first order term in the Taylor series. Thus, when £ < W we will have 



AQ 2 



2 i)2 



m z n £ 



+2V2^(Q 2 t ) 



5 2 n + 2W 2 ) A 



1-e 



-1 2 /2{S 2 Q +2W 2 



(21) 



where ttiq = C2- 



3. Results 



We use the previously published 350- fim SHARP polarization map of OMC-1 of lVaillancourt et al 



2008 to determine the turbulent correlation lengths of the Stokes Q and U parameters and the un- 
polarized emission N . Therefore, we refer to Equation (|2T1) derived for the structure functions 



(AK 2 (£) 



m 2 K f 



4 



(5 2 K + 2W 2 ) A 



-f 2 /2(S 2 K +2W 2 



(22) 



in which K stands for Q, U, or N. The coefficient m 2 K belongs to the first order term of the Taylor 
series, and 5k represents the turbulent correlation length scales of Q, U, andN. The calcu lations 
are done using the SHARP beam radius W = 4'.'7 (or FWHM = 11"; see lHoude et al.lbood ). 
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FiguresQ] [2 and [3] show the results of our fits to the data for Q, U, and N, respectively. In 
all three figures, the data are shown with symbols. In the top graphs, the broken curve does not 
contain the correlated turbulent term of the function 



2V2ir(K?) 
but only includes the following two 



S 3 K 



(5 2 K + 2W 2 ) A 



-£ 2 /2(5 2 K +2W 2 ) 



(23) 



m 2 K i 2 + 2 



2tt(K 2 ) 



4 



(d 2 K + 2W 2 ) A 



(24) 



We fitted Equation (|24p for values of 0'5 < £ < 0f8 (more details can be found in Appendix A of 
Houde et al.l (|2009l )).Ttie intercept of the fit at £ = shows the level of the turbulent component 
(i.e., the second term in the above equation) that can be compared to the first term, which depicts 
the contribution from the ordered part when £ < W. 

The middle plot shows the same information only plotted as a function of i. We subtract the 
data points from the fit of Equation (|24p and show the results in bottom graphs (symbols). The 
correlated turbulent term (Equation |23| ) is fitted to the data (broken curve) with 5 k as the only 
fitting parameter to match the width of the function. Even though the autocorrelation functions 
of Q, U , and N are assumed to have Gaussian patterns and the data points in the top and middle 
graphs in all three figures follow the fits quite well, it is quite unlikely that these are realistic models 
for these functions. The solid curve in the bottom graphs shows the contribution of the telescope 
beam alone (i.e., when 5k — in the exponent of Equation |23|). 

The results from the fits are tabulated in Table [TJ We have measured the turbulent correlation 
length scales of Stokes Q and U, and N to be approximately 9"4 ±0^1, 7"3 ± tf'l, 12"6 ± Of! 2 (or 
20.5 ± 0.2, 16.0 ± 0.2, and 27.5 ± 0.4 mpc at 450 pc, the adopted distance for OMC-1). 



Table 1. Results from our fit of Equation (|24|) to the dispersion data for OMC-1. 









9'.'4 ± Of! 1 


T!3 ± 071 


12 / .'6 ± 0".2 


20.5 ± 0.2 mpc 


16.0 ± 0.2 mpc 


27.5 ± 0.4 mpc 
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Discussion 



The study and analysis presented in this paper seek to characterize the magnetized turbulence 
in star forming regions in molecular clouds through the determination of the turbulent correlation 
length scale 5. 

Dust grains are present in a variety of astronomical environments, such as molecular clouds, 
and also tend to couple to the magnetic field threading these regions mainly through their magnetic 
moment. Comparison of the turbulent correlation length scales of different components of the 
dust emission with the results from other techniques and methods is helpful in achieving a better 
understanding of turbulence, magnetic fields, and their interactions. The turbulent correlation 
length scales are evaluated from the autocorrelation function of the emission, which is the Fourier 
transform of the emission power spectrum, and therefore 5 being inversely proportional to the width 
of the power spectrum, yields valuable information. The turbulent correlation length scales of the 
Stokes parameters of the polarized emission (5q ~ 21 mpc and Sjj — 16 mpc) a r e in a greement with 
that of the turbulent magnetic field (5 ~ 16 mpc) determined by lHoude et al.l (120091 ) . 



Another comparison can be made using the results presented in iLi &: Houdd (|2008n for the 
turbulent power spectrum of coexisting ion and neutral molecular species in M17. The spectra for 
these two species share the same pattern in the inertial range, which is expected due to the tight 
coupling between ions and neutrals through flux freezing. However, their spectra cease to follow a 
common pattern as they decouple through turbulent ambipolar diffusion . The values determ ined 
for 5q and 5jj in this paper are consistent with the analysis presented in iHoude et al.l (120091 ) . i.e., 
both 5q and 5u are larger tha n the ambipo l ar di ffusion scale, 5 ad, which was recently measured 
to be 9.9 mpc in Orion KL by IHoude et al.l (|201ll ). These measurements are also consistent with 
more general theoretical expectations related t o turbulent ambip olar diffusion (jLazarian et al.ll2004l : 



McKee k Ostriker I200I: iFalceta et al 
ments \\A & Houdel |200 



2010 



iHezareh et al 



2010|). 



Tillev et al.ll2010l ) and other observational measure- 



On the other hand, there is a measurable difference between the estimated value for 5m (— 28 
mpc) and the correlation length scales of the Stokes parameters, i.e., 5q and Sjj- This result implies 
that the power spectrum of N is different from that of Q and U. A lack or weakness of coupling 
between the magnetic field and the dust particles responsible for the unpolarized part of the emission 
can therefore be inferred from the results. 

The Stokes parameters Q and U are each derived from simultaneous measurements of two 
orthogonal polarization states. Considering the uncertainties in our estimates of 5q and Su, the 
difference between the values of these two paramet ers is significant and reveals the presence o f 



anisotropy, which is expected for bo th incompressible (Goldrei ch fc Sridharlll995l:ICho et al.l (|2002l )) 



and compressible MHD turbulence ( Cho fe Lazarian 2003 : Kowal fe Lazarian 201ol ). More precisely. 



it is predicted that the autocorrelation function of magnetized turbulence will have a longer length 
scale in a direction parallel to the magnetic field (as compared to an orientation perpendicular to 
the field). We have measured the mean polarization angle in OMC-1 to be approximately 30°, 



-9 - 



as can be visually verifi ed with the polarization map of this region presented in Figure 1(a) of 



Vaillancourt et al.l ((2008). This implies that the Stokes U parameter will be dominated by the 
emission polarized (approximately) along the mean polarization vector and therefore perpendicular 
to the mean magnetic field. We thus expect that it will have a shorter correlation length scale 
compared to the Stokes Q parameter, which will more or less be equally representative of emissions 
along and perpendicular to the magnetic field. This is consistent with our results, as presented in 
TableU 

Obtaining a complete turbulent power spectrum is only possible through high resolution ob- 
servations and sufficient sampling in space, which are not available to us through polarization 



measu rements with SHARP. Such analyses have, however, been recently conducted bv lHoude et al 



5. Conclusions 



In this paper we take advantage of the methods presented in Hildebrand et al. (j2009l ) and 



Houde et al.1 (|2009l ) to determine the turbulent structure functions of the Stokes parameters ( Q and 



U) of the polarized emission and unpolarized emission N. Subsequentl y, the solutions are fitted to 



the previously published 350-/zm SHARP polarization map of OMC-1 of lVaillancourt et al.l (|2008l ) to 
estimate the turbulent correlation lengths of the Stokes parameters 5q and 5jj and the unpolarized 
emission 5m- Our results are consistent with the results presented in Houde et al. (2009, 2011) and 
may indicate presence of anisotropy in the magnetized turbulence. We also infer a weak coupling 
between the dust component responsible for the unpolarized emission N and the magnetic field B 
from the significant difference between their turbulent correlation length scales. 
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ment Fund programs. The CSO is funded through NSF AST 05-40882. This work has also been 
supported in part by NSF grants AST 05-05230, AST 02-41356, and AST 05-05124. 
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Fig. 1. — The structure f unction ( AQ (^) ) usm g t ne previously published 350- /im SHARP polar- 
ization map of OMC-1 of IVaillancourt et al.l (120081 ) . Top: Equation ()24j ) (broken curve) is fitted to 
the data (symbols) and plotted as a function of £ 2 ; middle: similar to top but plotted as a function 
of £; bottom: the difference between the data points and the fit of Equation (j24"j) is shown as symbols 
with the broken curve depicting the fit of Equation (|23p with 5q = 9 / . / 4±0 / .'l. The solid curve shows 
the contribution of the telescope beam alone (i.e., when 5k— 0). The telescope beam is assumed to 
be Gaussian. 
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Fig. 2. — Same as Figure 1, but for ^ All 2 (£))■ The turbulent correlation length is measured to be 
5 u = 7"3±0'!l. 
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Fig. 3. — Same as Figure 1, but for (^AN 2 (£)^. The turbulent correlation length is measured to be 
5 N = 12% ± 0'/2. 



